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We study the effect of random forces on a double stranded DNA in unzipping the two strands, 
analogous to the problem of an adsorbed polymer under a random force. The ground state develops 
bubbles of various lengths as the random force fluctuation is increased. The unzipping phase diagram 
is shown to be drastically different from the pure case. 
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Biological processes such as DNA replication and RNA 
transcription get initiated and then proceed by unzip- 
ping of double stranded DNA (dsDNA) by various en- 
zymes like helicase, polymerase etc.Q. These enzymes 
exert force on dsDNA often directly, but also indirectly 
by maintaining a fixed distance between the DNA strands 
(a fixed distance ensemble). It was predicted theoret- 
ically that a dsDNA unzips to two single strands (ss) 
when the force exceeds a critical value which depends 
on temperature 0, d, 0, S, S 01- The unzipping transi- 
tion have also been studied experimentally by using single 
molecule manipulations Q. A consequence of this force 
induced unzipping transition is that in a fixed distance 
ensemble, a dsDNA shows a coexistence of a zipped (ds) 
and an unzipped phase, known as a Y-fork in biology, 
with a "domain wall" separating the two phases as the 
junction of the Y. The motion of the domain wall under 
a local instability either by a direct force or by the motor 
action of the enzyme leads to a gradual nonequilibrium 
unzipping (as a propagating front) Q. However for RNA 
polymerase or even dnaA helicase, there is an additional 
need to open up a bubble at the right place ("origin" 
for initiation). Though some force induced mechanism is 
expected here, but not much details are known. 

The unzipping of DNA is a competition between the 
binding of the base pairs (to be called monomers) and the 
orientation of individual links connecting the monomers. 
A force applied at any point (say the end point) gets 
transmitted to individual bonds to orient each one in the 
direction of the force. However, for a real DNA, there 
are various sources of inhomogeneities. Commonly stud- 
ied cases of sequence heterogeneity, stacking energy, etc. 
generally affect the bound or the ds part of the DNA. In 
a cellular medium, there arc single strand binding (SSB) 
proteins which bind to single strands. Such bound ob- 
jects can lead to variation in the response of individual 
links (or bonds) to the external force. Wc study the role 
of such binding proteins in DNA unzipping by consider- 
ing a situation where the nature of binding is modeled by 
a randomly oriented force. To avoid many independent 
randomness, we avoid other heterogeneities like sequence 
heterogeneity We compare our results with the DNA 
pulling at the end. 



Polymers with various types of randomness or disorder 
constitute a special class of disordered problems because 
of the occurrence of non-symmetry related (i.e., config- 
urationally distinct) degenerate ground states[ll|. The 
barriers (e.g., in space, in energy) separating these states, 
the widths of the local wells, etc., then determine the 
equilibrium and also the dynamic behavior of the poly- 
mer. The random problem we discuss here is not one of 
those studied earlier, but rather shows certain unique fea- 
tures, notably degeneracy of the ground state at special 
points. All these features, biological motivations apart, 
make this random force problem stand out on its own. 
A close cousin is the problem of a directed polymer in 
a random medium. But this problem is controlled by a 
"T = 0" fixed point (in a renormalization group sense) 
in Z? = 1 -|- 1 but a disorder induced phase transition oc- 
curs in higher dimensions, D > 3. In the present case we 
shall see a disorder-induced transition, even in D = 1 + 1 
dimension. 

The strands of the DNA are complementary to each 
other and every base in one strand knows its pair on the 
second strand. Considering the case where the random- 
ness works similarly on the two strands (i.e., either trying 
to keep the strands closer or unzig), the Hamiltonian in 
the continuum can be written ast 
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where ri{z) is the d-dimcnsional position vector of a 
monomer at a length z along the contour of the ith 
strand from one end z = 0, A'^ is the length of each 
strand or polymer, V{r) is the binding potential, r{z) = 
ri(z) —r2{z), and g(z) is a random force. Both the poly- 
mers are tied at end z = 0. For the "pure" problem, 
g(z) is constant and the force term reduces to the stan- 
dard form — g • r{N). The first two terms on the right 
hand side represent the elastic energy or the connectivity 
of each polymer (taken to be Gaussian). The base pair 
interaction is for monomers at the same location on the 
two strands. It follows that an equivalent description can 
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be obtained in which the two strands of the DNA are re- 
placed by a relative chain. The Hamiltonian in terms of 
r(z)is 
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with appropriate rescaling to make the elastic constant 
unity. This H also describes the problem of peeling of 
an adsorbed polymer by a pulling force 
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Natu- 
rally, both the problems have similar universal behaviour 
and many features of DNA unzipping (at least quali- 
tatively) can be understood by studying the unzipping 
of an adsorbed polymer. It is to be noted that if z is 
taken as an extra dimension (albeit different from the 
other d spatial coordinates), then these Hamiltonians 
also represent directed polymers in D = c? -I- 1 dimen- 
sions. This is the description we use in this paper. We 
choose uncorrelated randomness with [g(z)]dis = and 
[gi{z)gj{z')]dis = g'^SijS{z - z'). Here, [•••]dis denotes 
the quenched average over force realizations. We write 
g(z) = gCiz) where the C represents the random direction 
with unit variance. 

The aim of the present paper is to study the effect of 
a random force of zero mean on a bound or adsorbed 
polymer below its thermal unbinding or desorption tem- 
perature. If one is away from the DNA melting point or 
the thermal desorption transition, then the characteris- 
tics of the unzipping transition of the pure problem is not 
sensitive to the dimensionality of the system, as seen ex- 
plicitly in the exactly solvable directed polymer problem 
in different dimensions including D = 1 -|- 1 0, 01 • With 
that in mind, we use the discretized directed polymer 
model in Z) = 1 + 1 dimensions. 
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FIG. 1: (a) Schematic diagram of a directed polymer in 
D = 1 + 1 dimensions with an attractive wall at x — 0. The 
direction of the random force gC,{z) is shown by the arrows on 
each bond, (b) g (in units of e) vs T (in units of fcs/e) phase 
diagram. The points are for the randomly forced polymer and 
solid curve is for the pure case. The same convention is used 
in all the plots. 

The lattice version of Eq. [2]is a polymer in = 1 + 1 di- 
mensions, directed along the diagonal of a square lattice 
(Fig. [ija)). At X = 0, there is an attractive, impene- 
trable wall with binding energy — e (e > 0) which favors 
adsorption of the polymer on the wall. For DNA this im- 
penetrability implies mutually avoiding chains. One end 
of the polymer is always kept anchored on the wall while 



the other is left free. On each bond between the two con- 
secutive monomers, there is a random force g{z) = gCi^) 
(^(z) same for a layer) which is always perpendicular to 
the wall. The magnitude g related to the standard devia- 
tion of the force is kept fixed but the direction ({z) = ±1, 
i.e., either towards the wall or away from it, is chosen ran- 
domly with equal probability so that the average force, 
[.9C(^)]dis = 0. By averaging over the force configurations 
on the polymer (quenched averaging), we find that even 
in the absence of a fixed pulling force at the end, there is 
an unzipping transition if the variance of the force fluc- 
tuation exceeds a critical value. The force-temperature 
phase diagram shows an increase of the critical force with 
temperature (see Fig. [ijb)). The critical force starts de- 
creasing only near Tc and becomes zero at Tc- This is to 
be contrasted with the pure adsorption problem where 
the critical force decreases monotonically with the tem- 
perature. 

For every realization of the randomness, the partition 
function can be calculated exactly as 
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where W = [l + {e^'^ — l)6x,o\ ■ The superscript a in 
above expression denotes a particular realization, and 
Zn{x) is the partition function for a polymer with the 
nth (or the last) monomer at x. Physical quantities are 
to be averaged over the realizations. Quenched averag- 
ing is relevant here because the time scale of changes in 
the source of randomness is much slower compared to 
the thermalization of the polymer. One may note that 
an annealed averaging of Eq. [1] would yield an effective 
pure adsorption problem without a force though with a 
reduced elastic constant. The quenched averaging is dis- 
tinctly different. 

To monitor whether the polymer is zipped or unzipped, 
one needs the average distance of the last monomer from 
the wall ((•••) denoting the thermal averaging) [(x)]^;^ = 

[J2x^ ^l^^i^)/J2x^N'H^)]dis , and the isothermal ex- 
tensibility, being the response to the force, can be 
expressed in terms of position fluctuation of the end 
monomer 
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In the presence of a fixed applied pulling force at 
the end, the recursion relation can be solved exactly^, 
H, [13, In this case, the critical force is gc{T) = 

{kBT/2) In \e^'^ — l] with a temperature driven classical 
second order desorption transition at Tc = e/ln2 for 
5 = [l3|(see Fig. [Ub)). Henceforth fcs is absorbed 
in the definition of T {T = kBT). The average distance 
of the last monomer from the wall remains zero for any 
g < gc{T) but becomes proportional to A'' as soon as 
g > gc{T). This shows that the force induced unzipping 
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is a first order phase transition. There is no reentrance 
unlike the dsDNA problem because of absence of zero 
temperature entropy of the bound state. 

For the random problem, the ground state is the com- 
pletely bound polymer for zero force and new states 
can be obtained by flipping the adsorbed monomers. 
Fig. mja) shows the four possible force configurations of 
a monomer which is adsorbed on the wall. Let rii, 712, 
and 714 be the numbers of such configurations, then we 
have ni + + JT-s + = N /2, since the geometry of 
the model permits only N /2 monomers on the wall. For 
small 5, there is no gain in energy in flipping. Only ver- 
tex for which we can gain is vertex (Hi), if — e+2g > — 2g. 
Therefore, if g > e/4, one expects to see small bubbles. 
Below this force, the ground state is unique. One may 
note that the last term of Eq. [1] on integration by parts 
contributes a fixed force at the end plus a force gradi- 
ent which acts locally on the chain. For a negative force 
gradient, the strands minimize the free energy by max- 
imizing the distance between them. Thus, the bubble 
formation is not restricted to lattice models only and has 
wider validity. For g > e/4, after flipping all the type 
(Hi) vertices, the average energy is 

N 

Eq = -{ni+n2 + n4)e-2g{n2 + n3) = - — (3e + 45), (5) 

o 

taking all the four possible vertices to be equally proba- 
ble. Equating this with the energy of the unzipped state, 
—Ng, in which the polymer favors a configuration where 
each of the bond gets oriented in the local force direc- 
tion, we get the critical value of the force, gc ~ 3e/4. 
This analysis shows that there is a critical force fluctua- 
tion above which the polymer favours the unzipped state 
at T = 0. Similar arguments would indicate that larger 
bubbles are possible for g > e/2. These bubbles are dif- 
ferent from the bubble (eye) phase observed when the 
pulling force is applied in the interior of the DNA(i^. 
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FIG. 2: (a) The energy cost in flipping a monomer on the wall, 
(b) Ground state energy per monomer, Ea/N, as a function 
of force obtained for TV = 256 using exact transfer matrix at 
T = 0.001 (circles). Estimates from single flip Monte Garlo 
at T = for = 1024 (square) (e = 1). 

By using the exact transfer matrix for the recursion re- 
lation of Eq. ([3]) (and using log's to increase numerical ac- 
curacy), the ground state energy can be determined from 
the free energy at very low temperatures. By averaging 



over 10^ force realizations at T = 0.001, the estimated 
ground state energy per monomer for N = 256 is shown 
in Fig. [2{b) for various values of g. The same quantity is 
also calculated using the single flip monte carlo at T = 
(squares) for N = 1024. Data from both the methods 
match excellently and are in agreement with the predic- 
tion by the above analysis (Eq. (O) for e/4 < g < 3e/4. 
The data also show that for g = e, the free energy is 
—Ng, which is the energy of the unzipped state. 
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FIG. 3: (a) (a::) vs g, (b) extensibility for four different samples 
of length — 1024 and the average over 10^ such samples 
(thick solid lines). (c) Typical configurations for A'^ = 128 for 
three different g as indicated, (d) [{2;)]^^^ vs A*" for g = 2e. 
The solid line is the best fit to the data. All at T = 0.3 with 
e = 1. 

A surprising feature of the model is that the ground 
state configuration of the polymer depends on g with 
degeneracies appearing when g is an integral multiple of 
e/4. For example, for g = e/4, flipping of vertex (iii) does 
not cost any energy. On an average, there are N/8 such 
two fold degenerate vertices. Therefore, the entropy of 
the ground state for g = e/4 is 5*0(5 = e/4) = (TV/8) In 2. 
Similarly, if in a configuration, vertices (i) and (iii) are 
side by side, then flipping of these two vertices do not 
cost any energy if g = e/2. Further opening of bubbles 
would be possible by taking advantage of the forces on 
sites away from the wall. One can similarly argue for 
other values of g. There is a gradual increase of bubble 
sizes as g is increased beyond g = e/4. From the free 
energy we calculated the entropy at low temperatures. 
After averaging over 10^ realizations and then extrapo- 
lating to T = one gets the zero temperature or residual 
entropy which for g = e/4 agrees nicely with the entropy 
calculated above. 

In order to explore the effect of thermal fluctuations, 
we study the force-distance isotherms. The (a;) vs. g 
isotherms for four different samples of length N = 1024 
and also the average over 10^ samples are shown in 
Fig. Elja). The isotherm of an individual sample shows 
steps similar (but different in origin) to the steps seen in 
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the force-distance isotherms when an adsorbed polymer is 
subjected to a pulhng force in a random environment [isj . 
The steps for the random medium case appear due to the 
pinning of the polymer in the attractive pockets formed 
by the random distribution of energy on the lattice sites 
with the force attempting to depin from these pockets. 
In the present case, there is a mutual competition be- 
tween the set of bonds with the force towards the wall 
and the others with the force in opposite direction. The 
former set would like to opt for a state with monomers 
on the wall while the rest trying to unzip them. The 
steps in the isotherm lead everytime to a comb-like ex- 
tensibility as shown in Fig. [3ljb) . This indicates that the 
polymer responds in a "jerky" manner, by opening up 
local pockets of pinned regions. This is corroborated in 
Fig. ini^c) , where the locations of individual monomers for 
a chain of length = 128 are plotted for three different g 
chosen from different plateaus of the isotherm. It indeed 
shows adsorbed regions followed by unzipped regions and 
unzipping of adsorbed regions as g increases. 

Figure IHtd) shows [(x)]^,,, vs N for .g = 2e at T = 0.3 
(unzipped region). The solid curve which is the best fit to 
the data shows that [{x)]^^^ = aN" with a = 0.753±0.006 
and V = 0.505±0.001, i.e. the polymer stays at a distance 
of ^/N from the wall and the configuration is not like a 
directed polymer in a random medium for which v = 2/3. 
One may add that in the random medium problem the 
polymer gets fully stretched by the force in the unzipped 
state (i.e., depins from all the pockets). 

The phase diagram is obtained from the disorder av- 
eraged extensibility. Fig. [DJa) shows the plot of ex- 
tensibility [x]dis 9' the polymer of lengths N = 
256, 512, 1024 and 2048 at T = 0.3 averaged over IQ-' re- 
alizations. The growth in the peak with N indicates the 
possibility of a divergence in the thermodynamic limit 
{N oo). Its position can be located using the finite 
size scaling of the form [x]dis — N'^y{{g — gc)N'^), where 
gc is the force at which the discontinuity is located in 
the thermodynamic limit and d and cj) are the charac- 
teristic exponents. The data collapse obtained using the 
Bhattacharjee-Seno procedure [15[ is shown in fig. |D(b) 
which gives d = 0.58 ± 0.02, (j) = 0.22 ± 0.01 and g^ 
1.46 ± 0.04, with d/(j) w 2.6, so that [x]dis \g ~ gc. 
Similar collapse is also found for [{x)]dis (not shown). 
This indicates a continuous transition though the expo- 
nent is different of that of the heterogeneous DNA se- 
quence. We adopt this procedure at various temperatures 
to trace out the g — T phase diagram. For temperatures 
near Tc, lengths up to = 8192 are used. 

The phase diagram is shown in Fig[T]Jb). In contrast 
to the pure case, there is an increase in gdT) because of 
the loss of the residual zero temperature entropy of the 
bound state 

In conclusion, we have introduced and studied the 
problem of a randomly forced DNA or its equivalent ran- 
domly forced adsorbed polymer. The fluctuating force 
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FIG. 4: (a) Extensibility [xldis vs g at T = 0.3 with e = 1 for 
various chain lengths, (b) Data collapse of the extensibility. 



unzips the DNA by a gradual increase of bubble sizes. 
For some discrete values of the force (integral multiple 
of e/4 in our case), the ground state of the random force 
problem is degenerate and contains bubbles, which is in 
contrast to the pure problem where there is only one 
ground state. It suggests the possibility of opening up 
local bubbles in selective regions (proper pockets of force 
distribution) without unzipping the whole DNA. Only 
experiments can tell us if the initiation of DNA replica- 
tion or RNA transcription is through such a mechanism 
with the smaller molecules playing the role of random 
force agents. 

R.K. thanks SINP for hospitality. 
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